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Abstract 


5 331 ! 

In this article a theory of the general perturbations of a planetary 
system is developed. 

The perturbations of the position vectors of each planet are ex- 
panded into series arranged in powers and products of the masses m. , 
m 2 , . . . , m n of the planets constituting the system. Perturbations of 
different orders are obtained in the form of series containing the purely 
periodic, the secular and the mixed terms, in accordance with standard 
astronomical practice. The influence of the lower order perturbations 
on the higher order ones is determined. 

Typical differential equations are formed determining those per- 
turbations of the i th planet which are proportional to m k , m k m p , 
m k m p m s The right sides of these differential equations are 

obtained as the corresponding terms in the Maxwellian expansion of the 
gravitational forces in terms of multipoles. The momenta of these 
multipoles are the perturbations of all possible orders. 

The explicit calculation is carried out here for the perturbations of 
the first, second and third orders, and the procedure for determining 
the higher order perturbations is outlined. 

Decomposing each perturbation of any particular planet nq along 
the undisturbed position vector r. , along the undisturbed velocity v. , 
and along the unit vector R. normal to the undisturbed orbital plane, 
we deduce the differential equations to a form easily integrable by 
quadratures. After the integration it is more convenient in practical 
applications to replace this decomposition of perturbations by the de- 
composition along r. , R. x 7. and R. . The problem of the constants of 
integration is treated for the case of the elements osculating at the 
initial moment and for the case of the mean elements. The results 
given here extend and generalize the author's previous results on the 
case of the whole planetary system. The method suffers, however, 
from disadvantages common to all astronomical methods of the general 
planetary perturbations. It is not applicable to pair of planets if 
their orbits approach each other very closely. 
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Sir. = the perturbations in the position vector of the i th planet 
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ON THE GENERAL, PERTURBATIONS OF THE 


POSITION VECTORS OF A PLANETARY SYSTEM 
INTRODUCTION 

In this article a theory of general perturbations of a planetary sys- 
tem is developed. The perturbations in the position vector of each 
planet is developed into a series in powers and products of the disturb- 
ing masses and into a series containing the periodic, the secular, and 
the mixed terms with respect to time. 

Such a way of representing the integrals of the disturbed motion is 
in accordance with standard astronomical practice. From the pur ely math- 
ematical standpoint, this solution can be affected by all the difficulties 
associated with the near resonance conditions caused by the small 
divisors . 

We establish the differential equations for perturbations propor- 
tional to the powers and products of masses in a form integrable by 
quadratures. The explicit calculation is carried out through the per- 
turbations of the third order. In our planetary system it is rarely 
necessary to include the perturbations of the fourth and higher orders. 
However, an outline of the procedure for including the perturbations of 
even higher orders is indicated here. 

The problem of a direct determining of the general perturbations in 
the position vectors, including the effects of higher orders, became pos- 
sible only in recent years with the advent of electronic computers. By 
decomposing the perturbations along the directions of r. , v. and R. , 
[Musen and Carpenter, 1963] one can integrate the variational equation 
of the problem by HilFs [1874] procedure directly without resorting to 
the method of variation of astronomical constants. We shall use here 
the same decomposition as an intermediary step; but the final decomposi- 
tion of the perturbations will be along r . , R. xr. and R i , in order to 
reduce the components of the disturbing term on the right side of the 
variational equation to a simple form. In computing the higher order 


I 


perturbations, it will be necessary to expand the disturbing forces in 
powers of perturbations in the position vectors. Maxwell’s method of 
expanding the electrostatic potential in terms of multipoles by employing 
symbolic operators [1904, 3rd ed] can be used profitably also in planetary 
theories. In our exposition the moments of the multipoles are the per- 
turbations of different orders of the position vectors. Evidently any other 
way of expanding the disturbing forces in terms of the perturbations of 
the position vectors will lead to a duplication of Maxwell's expansion, 
but through a more laborious writing. 

O O 


In the theory of perturbations of the position vectors we achieve 
economy of theoretical thinking as well as economy of computing 
machine time, because a set of homogeneous operations is being con- 
stantly repeated. All these circumstances suggest that future methods 
of calculating general perturbations will be based on the expansion of 
the perturbations in the position vectors directly. 


The Differential Equations of the Problems 

Putting 


p - r - a 


we shall make use of Maxwell expansion of the spherical functions as 
defined in terms of multipoles. We have 


<£( n ) 




r 

j l.3---(2n-l) 

[ p 2n + l 


n 



k = l 


k = 1,2, • 

n = 1,2,3 


n 


1.3 •• • (2n - 3) 
p 2 "-i 


1.3 ••• (2n - 5) 

p 2 "" 3 


L 

L 


n 



k = 5 


(i) 


2 


♦ 



here a . a^ , * . . , a are constant vectors, V is the del-operator with 

1 2 n 

respect to r, and the sums 


n 



I I 
k= 3 


designate the sums of all terms as obtained from the first term by 
means of the permutations of all n indices. In particular, we have 
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The gradient of the spherical function <p( n > is obtained from the 
Maxwellian expansion of ^ n + i ^ simply by replacing the moment a n+1 
by the idemfactor I. Thus from (Z) - (6): 
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The differential equation of the disturbed motion of the i th 
can be written in the form 
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Taking the equation 
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into account and introducing the vectorial differential operators 


E>. = V. exp(5r. • V. ) , (12) 

D i . = V. exp(8p..-V.) , (13) 


which perform the development of (11) into a power series in the com- 
ponents of S? k (k = l f 2, ••• n) , we obtain from (11) the differential equa- 
tion for Sr. in the form 



f: 




(14) 


The perturbation vector 57. can be developed into a series with 
respect to the powers and products of the disturbing masses. We put 


s7. = — y r a + — T tV + 1 y T a py + ■ 
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where 7^ is proportional to m a , and 7^ is proportional to m a , etc. 
The factors in front of the sums in (15) are introduced in order to re- 
move large coefficients in higher approximations. 


We define 7*^ , 7^^ , . . . in such a way that they remain invariant 
under the permutations of the upper indices: 
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etc. 


Symbolically, 


where in performing the development and the symbolic '’multiplications ' 1 
the indices are not being added, but written in a row. We have also 


r. - rf 


0 . 


We shall now deduce the differential equation for determining per- 
turbations of the form 


■Tk -pep fkpq rkpqs ... 


(i, k, p, * • ■ = 1, 2, •• • n) , 


first under the assumption that there are no identical indices among 
k, p, q, . . . Retaining only the substantial terms, we have 
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Substituting these values into 
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and again retaining only the necessary terms, we deduce: 


D. - V. = ? k -V. V. + ( 7 kp - V. V. + 7 k -y. r ?• V. V.) 

i i iii v i ii a i i i i y 


( 18 ) 


+ ( 7 kpq • V. V. + 7 kp - V. 7 q -V V + r pq -V. r k • V V 

v i ii i iiii i iiii 


7 qk -y. 7 p -V. V. + 7 k - V 7 p -7. 7 q -7. V i 

i iiii liiiiii' 


+ |7 kpqs • V : V. + 7 kpq -V ; 7*-V ; V. 


11 1 1111 


+ 7 kps -v. 7 q -y. y. 4. 7 kq s -v 7 p -y y 

i i iii- i iiii 


7 pq S .y 7 k.y_ v. + (7 kp -y 7 qs -y y 

i iiii v i ii ii 


+ 7 kq -v 7 ps -y. y + 7 ks -y 7 pq -y y.) 

1 11 11 1 11 i x y 


+ — (7 kp -y 7 q -y 7 s -y v A r k<i -y r s -y7 p -y y 

2 v i iiiiii’i ii iii 


+ 7 ks -y 7 p -y 7 q -y y +7«-y 7 k -y. 7 s -y. v. 

l i i i iii i iiiiii 


7 ps -v 7 k -y 7 q -y v + 7 qs .y. 7 k -y. 7 p -y. y.> 

i 111111 i i i i i i i ' 


+ r k -V. "r p *V. ~r^ -V. 7 S -V. V. 

ii ii ii iii 


8 



(19) 


D. . =7. 

i i j 


+ P 


j i j 


7 ; - (p kp -V. 


V. + 
) 


Jv 

^ . . 
J J 


V 


(Tkpq ,y v. 

v ‘ J 1 J J 


P 


kp *v. 3s . *7. v. 


;pq.y. 

j i j 


•V. 7. 

j j 


i ^ 


-V. V. 

j j 


+ p 


i i 


•7. 


p p -‘ 

j i 




•y. v.) 

j 3 7 


These developments of the operators D - V. and D. . permit one to 
compute the general planetary perturbations even up to the fourth 
order if necessary. 

We will establish the differential equations for determining the 
perturbations up to the third order. In our solar system occasions re- 
quiring the perturbations of the fourth order will be, probably, very 
rare. However, in some cases of very sharp commensurabilities of 
mean motions, the question remains open and further numerical in- 
vestigations are necessary. Substituting (18) and (19) into (14) and re- 
taining only the typical operators, we have: 
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Comparing the terms of the same degree in the disturbing masses in 
the left and right sides of (20), we obtain the basic equations for de- 
termining the general perturbations up to the third order: 
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The equations (Zl) - (Z3) can be written in the form 
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The terms in the right sides of (27) - (29) are the partial gradients of 
the sums composed of the elementary spherical functions, with the 
moments equal to the perturbations in r. and p (i , k = 1 , 2, ... ,n) . 
Making use of (24) - (26) we obtain 
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also the expanded forms of the typical terms in F* p : 
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and the expanded forms of the typical terms in F kpq : 
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In a similar way the expressions for F^ pqs and for its typical terms 
can be formed* The process can be continued as far as necessary, 
until the perturbations become negligible* The restriction that indices 
k, p, q, . . . be different now can be removed. 


Integration Procedure 

The typical differential equation which appears in the problem of 
the general perturbations in the position vectors has the form 
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ana F* is a series with the periodic, secular and mixed terms. Equa- 
tion (40) is the variational equation of the two body problems. 

In order to reduce the solutions of (40) to quadratures we shall 
make use of the substitution 


(s. + 2w.) 7. - v. (2s, T 3 w . ) dt 

v i i J i i v i i / 
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This differs from the substitutions used by the author [1963] previously 
The substitution (41) was chosen because of a simpler form to which 
(40) is reduced thereby as compared with the earlier exposition. It is 
unnecessary to retain the index i in the further exposition: we can 

now omit it without loss of clarity. 

It follows from (41) that 
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Substituting (41) and (43) into (40), we obtain the vectorial differential 
equation 
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which can be integrated by quadratures. Forming the dot products of 
(44) and 
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We have from (46) 


w = K 3 + B, 


where we out 


(48) 


B= f— na - F'Rxr dt ; (49) 

j 

K 3 is the additive constant of integration. The integration of series in 
(49) is performed in a formal manner. 

The equation (45) can now be integrated using Hill's procedure 
LI 874 j . We obtain: 
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The vector f is considered as a temporary constant and it is replaced 
by r after the integration is completed. The integrand is a trigono- 
metrical series in the mean anomalies of planets and the auxiliary 
mean anomaly *£ associated with r , and it can also contain the purely 
secular and the mixed terms. After the integration is performed 4 
has to be replaced by l . Integrating the second integral in (51) by 
parts and replacing rxT by zero when it appears outside the integral 
sign, we obtain 

fib' R.fxf dt = - I ihR.vxt dt (52) 

J dt 2 j Gt 
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From (46), (b I ) and (52) we deduce: 
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As in (49), the integration is performed in a formal manner. 
In a similar way we obtain from (47) 
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where K 4 is the additive constant, of integration. 

The forms (30) - (39) of the disturbing terms in the right sides of 
the variational equations require the decomposition of F* , r K? , 
kpq , ... in the moving frame R.xF, R rather than in the frame 
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The decomposition (57) was suggested by Pop o vie [I960, 1961] and 
independently by the author ll 964j . In Popovic ’ s work, the perturba- 
tions of the first and of the second order in <f, r, £ are determ : d. 

The disturbing vector F can also be decomposed along r, Rx r , and 
R. Perhaps from the computational standpoint this decomposition is 
the simplest one; it appears in several theories of the general planetary 
perturbations , either directly or as an intermediary step. Let us put 
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and replacing v by its decomposition as given above we have, after 
some easy vectorial transformations. 
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M is a sine series in Z and b , and N is a cosine series in the same 
arguments. In order to obtain these series we have to obtain the 
development of 
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The equation (49) can be written as 
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The computation oi integrands in (4 
development oi some other express 
of 1 /p..; of powers of s k /r k ; oi the 
and the triple oroducts R. * r.x r . . 
easiest way to obtain all these dove 


9 5 ), ( 53 5 ) and (55) requires also the 
ions: for example, of the odd powers 

scalar products f. * r. and r * R . ; 
It scorns that the simplest and 
1 contents is by means of the single 


and double harmonic analyses. The lormulas 


r . ~ A . £ cos € . - e . j -r B i s in € . , 


e . s 


, rr . 772 


serve as a start. The representation of /-v in terms of the vectorial 


elements is very useiul it 


1 performing the double harmonic analysis. 


Determination of the Constants of Integration 


Each perturbation 
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itro 


duces six constants of integration. We shall designate them by 


K k . , K k p , K k ? q , . . 
j 1 3 * 3 1 


i,k,p,q, ... “ 1 , 2 , . . . , n 


j =1,2, . . . 6 


correspondingly. Consequently, the problem ot determining the gen- 
eral perturbations by developing them into power series with respect 
to the masses introduces an ini in it e set of constants ot integration. Of 
course, with the increasing number of the upper indices these con- 
stants decrease rapidly. 


In the final exoressions for the 


irbations, the constants are 


combined to form a set 01 only on independent constants; out mis is not 
seen explicitelv in a numerical planetary theory. At eacii step me con- 
stants must be determined separately from some additional conditions 
imposed on the elements or ircm tne initial position and velocity vectors. 


the mean elements 

aments osculating at the initial moment of time. In the case 


Two types of elements are being commonly used: 


ana tne 

of the mean elements the perturbations ci the true longitude in the 
orbital plane, as defined by these elements , shall not contain tne constant 
term, the purely secular term and the terms with periods equal to tne 
period of revolution 01 the planet. The perturbations m tne "third 

normal to the orbital plane, shaii not co mam the terms with 


coordinate , 


periods equal to the revolution 01 the planet. 
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e I e m e n t s , h ow ever, 
the eccentric, true, 
to he used in develop 


are not unique. They depend upon the choice of 
or mean anomaly as the basic independent variable 
ing the perturbations of a given planet. 


Hansen L1857J made use of the eccentric anomaly in his theory of 
minor planets* Hill Li 87 4 j preferred the true anomaly. In both 
choices we gain speed of convergence of series giving the perturbations 
of the first order, in comparison to the choice of the mean anomaly. 
However, the road to computation of the higher order perturbations of 
the whole planetary system by either of these two choices will be 
blocked so effectively that all the gains in the first approximation appear 
to be negated by the difficulties encountered in computing these higher 
order perturbations. For this reason the use of the undisturbed mean 


anomalies 


and , c onsequently , of th e unive r s al var table . 


time, is highly recommended in the planetary theories. This is done 
already by Hansen and by Kill in their theories of Jupiter and Saturn. 
In connection with this statement we say that the elements are mean it 
there are nc terms of the form 


K Q1 Kjt, K< c) cos b, Kb s > sin b 

in the perturbations of the true longitudes with respect to the undis- 
turbed orbit planes, and there are no terms of the form 

cos i s ) sin i- 


in the : ! third coordinates 1 ' H 

We express the perturbations in the true longitude in terms of the 
perturbations along r°, Rk?, and R. Taking into account the equations 
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We can write 
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these iwo last equations with respect 
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Substituting 
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we obtain into (63) 
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Putting: 
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And substituting these values into (64) - (66) and taking (59) into con- 
sideration we obtain 


X 

n . V 1 -e 2 


3 / 3 i 


n. \ r : 


(- K*. + e* K*.) rtj t 


(70) 


K* 


2 sin f . + — e . sin 2 f . 

^7772 V 1 2 1 


Kj . Vl-e 2 / , , 

+ — I 2 cos f . + — e . + — e . cos 2 f . 

n \ 1 2 1 2 1 1 




k* p o / a . \ 2 

1 = — (— ) (- K 1 *? + e. K!f?) n. t 

n. VTT^ n i V r i / 1 31 1 lW 1 


- 2 sin f. + — e. sin 2f 

n. V 1 2 * 


■) 


K*(P V 1-e 2 


n. 


*2 i ' /of 1 1 

— 2 cos f .. + — e . + — e . cos 2 f 


■) 


a 2 


+ -fi? + w*? 

r 2 4i i 


(71) 


27 


A* 


kpq 


n. YT^T n i V r « 


a . 


(- K ‘,r + *ir) n . t 


Vkpq 
1 i 


1 


, 2 sin f. + — e. sin 2f. 

n. rr^T v 1 2 - 


K*?* Vl-e 2 


2 i 


rl-e l 11 \ 

L ( 2 C ° S f i + 2 G i + 2 G i C ° S 2f V 


a 2 


+ — + 1^ 
_ 2 4 1 i 


r 2 

1 


(72) 


We shall make use of the developments 
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The coefficients in these developments are computed either by some 
analytical classic procedure or by means of a harmonic analysis if the 
eccentricity is not too small. 

The terms of the form 


K n ., K. t, K( c) cos 'i, K( s > sin Z 


(73) 


must be absent in the developments of (70) - (71). We have, keeping 
only the substantial terms. 


w i = a oi + ^Oi n i 4 + a u cos + /fin sin l i -r . . . , 
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In the process of computing the coefficients a and /3 , the machine re- 
jects automatically all the useless terms, unless we decide to obtain a 
complete development of the perturbations in the true longitude. The 
conditions for the absence of the terms (73) in the right sides of (70) - 
(71) lead to the equations: 
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Separating in Z k , Z kp , Z kpq , . . . the terms with the argument we have 

“ c ii cos ^ + s ii s * n ^ + • • • > 

Z kp - c k ? cos 'i T s kp . sin Z + . * . , 

ill 1 1 ’ 
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Z k pq - c k pq cos J r s ^ pq sin 't + . . 

i 1 1 1 1 


(i * k, p, q, ... =1,2, ... , n) 


and the conditions for the absence of terms with the argument 'i in 
£ k , £ kp , 4 kpq , . . . lead to 


K s c i 1,X + c x = 0 f (81) 

K 6 S*’ 1 + Sl = 0. (82) 


The lower index i and the upper indices k, p, q, , . . are omitted in 
(77) - (82). Evidently, this is not causing any ambiguity. 


Conclusion 

The results given in this article represent the extension and com- 
pletion of the results given in the author's previous articles on this 
subject. The theory given here can also be considered as a modifica- 
tion and generalization of Hill's planetary theory, with the latter's 
inconvenience s removed. The interdependent constants of integration 
peculiar to Hill's theory do not appear in the present exposition. The 
solution is given in the form which permits us to write immediately 
the differential equation for the general perturbations proportional to 
any prescribed product of masses. Moreover, the vectorial formalism 
permits us to penetrate into the structure of higher orders effects with- 
out great difficulty. Programming is also facilitated by the repetition 
of the homogeneous operations. 

The formula (57) permits one to obtain easily the decomposition of 
r i> r ] ° > • . • • along the axes of the inertial systems, if it is considered 
as necessary. 

On the basis of experience obtained at Goddard Space Flight Center 
one might expect that computing the perturbations of a given order for 
one planet will require only a few minutes. Considering the simplicity 
of the methods for the general perturbations in the position vectors, one 
is inclined to believe that such methods will consitiute one of the prin- 
cipal approaches to the problem in the not too distant future. 
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